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ABSTRACT 

Context. Accurate and repeated photometric follow-up observations of planetary-transit events are important to precisely 
characterize the physical properties of exoplanets. A good knowledge of the main characteristics of the exoplanets is 
fundamental to trace their origin and evolution. Multi-band photometric observations play an important role in this 
process. 

Aims. By using new photometric data, we computed precise estimates of the physical properties of two transiting 
planetary systems at equilibrium temperature of ~ 2000 K. 

Methods. We present new broad-band, multi-colour, photometric observations obtained using three small class telescopes 
and the telescope-defocussing technique. In particular we obtained 11 and 10 light curves covering 8 and 7 transits of 
HAT-P-23 and WASP-48 respectively. For each of the two targets, one transit event was simultaneously observed 
through four optical filters. One transit of WASP-48 b was monitored with two telescopes from the same observatory. 
The physical parameters of the systems were obtained by fitting the transit light curves with jktebop and from 
published spectroscopic measurements. 

Results. We have revised the physical parameters of the two planetary systems, finding a smaller radius for both HAT- 
P-23 b and WASP-48b, Rb — 1.224 d= 0.037Rj up and Rb = 1.396 ± 0.051 Rj up , respectively, than those measured in 
the discovery papers (Rb = 1.368 d= 0.090 Rj up and Rb = 1-67 d= 0.10 Rj up ). The density of the two planets are higher 
than those previously published (pb ~ IT and ~ 0.3 pj up for HAT-P-23 and WASP-48 respectively) hence the two Hot 
Jupiters are no longer located in a parameter space region of highly inflated planets. An analysis of the variation of the 
planet’s measured radius as a function of optical wavelength reveals flat transmission spectra within the experimental 
uncertainties. We also confirm the presence of the eclipsing contact binary NSVS-3071474 in the same held of view of 
WASP-48, for which we refine the value of the period to be 0.459 d. 

Key words, stars: planetary systems - stars: fundamental parameters - stars: individual: HAT-P-23 - stars: individual: 
WASP-48 - techniques: photometric 


1. Introduction 


Among the almost 2000 extrasolar planets known to exist, 
those that transit their parent stars are of particular in¬ 
terest. In contrast to the exoplanets detected with other 
techniques, most of the physical and orbital parameters 
of transiting extrasolar planet (TEP) systems are mea¬ 
surable with precisions of a few percent usin g sta ndard 
astronomical methods fe.glSeager fc MallemO rnelas 2003; 
Sozzetti et al.l 2007: ISouthworth et al. l2007 ~ lTorres et ahl 
20081: ISouthworthll2008[ ). Obtaining estimations of the plan¬ 
ets’ masses and sizes can give some hint and direction in 
discriminating between gaseous and rocky structure and 
therefore infer their formation and evolution history. In par¬ 
ticular, precise measurements of planetary sizes can pro¬ 
vide strong constraints for those theoretical models that 
try to explain the inflation mechanisms for highly irradi¬ 


ated gaseous planets. Indeed, after the discovery of a group 
of inflated planets (e.g. WASP-17b, TrES-4b and HAT-P- 

32. 1 Anderson et al1l2Ql it ISozzetti et al 120151: ISeeliger et al.l 

l2014f h several theories invoking, tidal friction, enhanced at¬ 
mospheric opacities, turbulent mixing, ohmic dissipation, 
windshocks or more exotic mechanisms, have been proposed 
to account for the slow cooling rate of these p lanets , resu lt¬ 
ing in a so unexpected large radius (see e.g. iBaraffe et al.l 
l2Q14l : ISpiegel fc Burrows Il20131: [Ginzburg fe Sari Il2015l and 
references therein). A deduction of their chemical composi¬ 
tion is also possible by looking for elemental and molecular 
signatur es in transmission spec t ra observed during transit 
events (|Seager fc Sasselovll2Q00l: lBrownll2Q0lf ). 

In this context, we are carrying out a large program 
to study these TEP systems and robustly determine their 
physical properties via photometric monitoring of transit 
events in different passbands. We are utilising an array 
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the projected spin-orbit angle through observation of the 
Rossiter-McLaughlin (R-M) effect. Their finding of A = 
15° d= 22° suggests an aligned and prograde planetary or¬ 
bit. A reanalysis o f the main paramet e rs of the system 
was perfo rmed by iRamon-Fox fc Sa~dal ( 2013f ). Recently, 
lO’Rourke et al.l (|2Q14l) reported accurate photometry of 
planet occupations observed at 3.6 /im and 4.5 /im with the 
Spitzer space telescope and at H and K$ bands with the 
Hale Telescope. They found the emission spectrum to be 
consistent with a planetary atmosphere having a low effi¬ 
ciency of energy transport from its day-side to night-side, 
no thermal inversion and a lack of strongly- absor bing sub¬ 
stances, similar to the case of WASP-19 b (|Mancini et al.l 
12013c). 


Fig. 1 . 3D plot of the known transiting exoplanets. The 
quantities on the axes are the planetary radius and the 
temperature, and their semi-major axes. The positions of 
HAT-P-23 b and WASP-48 b are highlighted by a blue and 
red point, respectively. Data taken from the TEPcat cata¬ 
logue. 


of 1-2 meter class telescopes, located in both of Earth’s 
hemispheres, and single or multi-channel imaging instru¬ 
ments. In some cases, a two-site observational strategy was 
adopted for simultaneous follow-up of transit events (e.g. 
ICiceri et al.ll2013l: iMancini et al.ll2Q13al l2014al) . We have so 
far r efined the measured parameters of several TEP system s 
(e.g. ISonthworth et al.ll2Qllt IMancini et al.ll2013al 120143). 
studie d starspot crossing events ([Mancini et al I 12013 d 
l2014b[ ). and probed opacity-induced vari ations of measured 
planetary radius with wave le ngth (e.g. Southworth et al.l 
l2012bl: IMancini et al.ll2Q13blR l2Q14bl : iNikolov et al.ll2Q13f) . 


In this work, we focus our attention on the HAT-P- 
23 and WASP-48 systems, both hosting a star with effec¬ 
tive temperature T e ^ ~ 6000 K and a close-in gas-giant 
transiting planet with a high equilibrium temperature of 
T eq ^ 2000 K. We have reported these characteristics in 
Fig.[TJ together with those of the other known TEPs (data 
taken from TEPCatQ). The other main properties of the 
two TEP systems are summarized in the next two sub¬ 
sections. Observations of 15 new transit events of the two 
planets, performed at two different observatories, are pre¬ 
sented in Sect. [2] together with the reduction of the corre¬ 
sponding photometric data. The analysis of the light curves, 
described in Sect.[3j is followed by the revision of the physi¬ 
cal parameters of the two planetary systems in Sect.0J We 
also present new light curves of the eclipsing binary NS VS 
3071474, which is located close to WASP-48. Our conclu¬ 
sions are summarized in Sect.[5l 


1.1. HAT-P-23 

The HAT-P-23 system, discovered bv iBakos et al.l (|201ll) . 
is composed of a GO dwarf star (mass 1.13 Mq, radius 
1.20 Rq and metallicity +0.13) and a hot Jupiter (mass 
2.09 Mj up and radius 1.37Rj up ), revolving around its par¬ 
ent star on a circular orbit with a peri od of P = 1.2 d. Using 
the SOPHIE spectrograph. iMoutou et al.l (|201lh measured 

1 TEPcat is the catalo gue of the physica l properties of tran¬ 
siting planetary systems (lSouthworthll20lIh and is available at 
http://www.astro.keele.ac.uk/jkt/tepcat/, 


1.2. WAS P-48 

WASP-48 is a TEP system composed of a slightly evolved 
F-type star (mass 1.19 Mq, radius 1.75 Rq and metallicity 
—0.12) and an inflated hot Jupiter (mass 0.98 radius 

1.67 Rj up and an orbital period of P = 2.1 d) ([Enoch et al. I 
l201lh . The lithium abundance and the absence of Call H 
and K emission suggests that the host star is old and evolv¬ 
ing off the main sequence. However, by measuring the age 
from th e rotat ion period, the star seems to be much younger 
([Enoch et al. I [2QlTh . This discrepancy can be explained if 
tidal forces between the planet and the parent star have 
spun up the latter, making it difficult to obtain a reason¬ 
able value for the age based on gyrochronology (e.g. iPontl 
2009]). The emission spectrum of W ASP-48 b was also inves¬ 
tigated bv lO’Ronrke et all ((2014) through IR observations 
of the planet occultations with the Spitzer and the Hale 
telescopes. The nature of the atmosphere of WASP-48 b was 
found to be very similar (moderate energy recirculation, no 
temperature inversion, absence of strong absorbers) to that 
of HAT-P-23 b. 


2. Observations and data reduction 

In this section we present photometric observations of eight 
transits of HAT-P-23 b and seven of WASP-48 b. For both 
systems, one transit was observed with a multi-band imag¬ 
ing camera. One transit of WASP-48 was simultaneously 
monitored with two telescopes from the same observatory. 
The details of the observations are summarized in Table[2j 

2.1. Calar Alto 1.23-m Telescope 

Seven transits of HAT-P-23 b and four of WASP-48 b were 
remotely observed using the Zeiss 1.23-m telescope at the 
Calar Alto Observatory in Spain. The telescope has an 
equatorial mount and is equipped with a DLR-MKIII ca¬ 
mera positioned at its Cassegrain focus. The CCD, which 
was used unbinned, has 4096 x 4096 pixels and a field of 
view (FOV) of 21.5' x 21.5' leading to a resolution of 0.32" 
pixel -1 . The telescope was autoguided and defocussed for 
all science observations. This observing mode consists of 
using the telescope out of focus to spread the light of the 
stars in the FOV on many more pixels of the CCD than 
normal in-focus observations. In this way it is possible to 
use longer exposures, greatly increasing the signal to noise 
ratio (S/N) and reducing the uncertain ties due to many 
sources of noise (Sou thworth et al1l2009l) . In our cases, the 
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Table 1 . Details of the observations presented in this work. 7V 0 bs is the number of observations, T exp is the exposure 
time, T 0 bs is the mean observational cadence, and ‘Moon ilium.’ is the fractional illumination of the Moon. 


Telescope 

Date of 

Start time 

End time 

iV obs 

T 

-L exp 

Fobs 

Filter 

Airmass 

Moon 

Aperture 

Scatter 


first obs 

(UT) 

(UT) 


( S ) 

(s) 




ilium. 

radii (px) 

(mmag) 

HAT-P-23: 
CA 1.23-m 

2011 08 19 

23:16 

00:52 

60 

80 

103 

Cousins R 

1.78 -A 

1.34 

67% 

13,70,90 

1.41 

CA 1.23-m 

2012 07 21 

00:34 

04:30 

139 

50 to 120 

101 

Cousins R 

1.07 -A 

1.78 

9% 

14,70,90 

0.69 

CA 1.23-m 

2012 07 26 

20:25 

02:37 

286 

30 to 130 

91 

Cousins R 

1.70 -A 

1.07 

58% 

23,70,90 

1.22 

C A 1.23-m 

2013 06 16 

22:50 

03:49 

120 

135 

151 

Cousins R 

1.81 -A 

1.07 

53% 

23,60,90 

0.95 

CA 1.23-m 

2013 07 03 

21:47 

03:49 

194 

85 to 100 

110 

Cousins R 

1.78 -A 

1.07 

15% 

16,70,90 

0.77 

C A 1.23-m 

2013 07 26 

20:25 

02:37 

87 

30 to 130 

61 

Cousins R 

1.70 -A 

1.07 

76% 

23,70,90 

1.59 

CA 1.23-m 

2013 07 31 

19:55 

03:02 

147 

145 to 155 

166 

Cousins R 

1.80 -A 

1.07 

28% 

25,70,90 

1.19 

CA 2.2-m 

2013 09 03 

19:39 

23:45 

113 

60 to 120 

133 

Thuan-Gunn u 

1.21 -A 

1.07 

2% 

10,25,40 

3.55 

CA 2.2-m 

2013 09 03 

19:39 

23:45 

117 

60 to 120 

133 

Thuan-Gunn g 

1.21 -A 

1.07 

2% 

23,55,70 

0.92 

CA 2.2-m 

2013 09 03 

19:39 

23:45 

118 

60 to 120 

133 

Thuan-Gunn r 

1.21 -A 

1.07 

2% 

21,45,65 

0.71 

CA 2.2-m 

2013 09 03 

19:39 

23:45 

116 

60 to 120 

133 

Thuan-Gunn z 

1.21 -A 

1.07 

2% 

17,50,70 

1.49 

WASP-48: 

Cassini 

2011 05 23 

22:04 

01:02 

91 

52 to 90 

92 

Gunn r 

1.30 -A 

1.04 

57% 

13,70,90 

0.82 

Cassini 

2011 05 25 

21:33 

02:45 

222 

50 to 90 

81 

Gunn r 

1.37 -A 

1.02 

38% 

13,70,90 

0.81 

CA 1.23-m 

2011 08 23 

22:56 

03:53 

259 

25 to 42 

63 

Cousins R 

1.77 -A 

1.27 

29% 

11,30,45 

1.90 

CA 1.23-m 

2012 09 12 

19:08 

01:25 

136 

120 to 160 

151 

Cousins / 

1.07 -A 

1.74 

5% 

20,70,90 

1.01 

CA 1.23-m 

2013 07 24 

20:10 

04:12 

214 

110 to 120 

125 

Cousins I 

1.05 -A 

1.59 

92% 

22,70,90 

0.61 

CA 2.2-m 

2011 08 23 

21:23 

04:40 

256 

50 to 80 

101 

Stromgren u 

1.07 -A 

2.36 

29% 

10,16,30 

2.88 

CA 2.2-m 

2011 08 23 

21:23 

04:40 

246 

50 to 80 

101 

Thuan-Gunn g 

1.07 -A 

2.36 

29% 

25,40,80 

1.19 

CA 2.2-m 

2011 08 23 

21:23 

04:40 

248 

50 to 80 

101 

Thuan-Gunn r 

1.07 -A 

2.36 

29% 

23,33,60 

1.09 

CA 2.2-m 

2011 08 23 

21:23 

04:40 

134 

50 to 80 

101 

Cousins I 

1.07 -A 

2.36 

29% 

18,28,50 

1.25 

CA 1.23-m 

2014 06 02 

20:34 

03:40 

184 

115 to 134 

126 

Cousins R 

2.12 -A 

1.05 

25% 

20,80,100 

0.89 


exposure times were fine-tuned at the beginning of each 
observation, together with the amount of defocussing, in 
order to properly optimize the S/N and have a maximum 
count per pixel for the target star between 25000 to 35000 
ADUs. Once the defocussing amount was set, it was kept 
fixed for the entire monitoring of the transit (a typical PSF 
of the target covered a region with a diameter from 15 to 25 
pixel). In some cases, it was necessary to modify the expo¬ 
sure time during the night to avoid the CCD saturation or 
to account for changes in counts caused by variation of the 
air mass of the target or weather conditions (e.g, dramatic 
variations of the external temperature or sudden appear¬ 
ance of cirrus and veils). The filter used to observe HAT- 
P-23b was Cousins R‘ for WASP-48 b the first and last 
transits were observed through Cousins R and the other 
two transits through Cousins I (Tabled . 

2.2. Cassini 1.52-m Telescope 

Two transit light curves of WASP-48 were obtained on 2011 
May 23 and 25 at the Astronomical Observatory of Bologna 
in Loiano, Italy. Observations were carried out with BFOSC 
(Bologna Faint Object Spectrograph & Camera) mounted 
at the Cassegrain focus of the 1.52-m Cassini telescope (see 
Tabled . The 1300 x 1340 pixels CCD has a FOV of 13' x 
12.6' resulting in a resolution of 0.58" pixel -1 . For both 
transits a Gunn-r filter was used and the telescope was 
autoguided and defocussed. During the acquisition of the 
science images, the CCD was windowed to decrease the 
readout time, thus achieving a higher temporal cadence. 

2.3. Calar Alto 2.2-m Telescope 

One transit of each object was observed with the Bo nn 
University Simultaneous Camera (BUSCA, iReif et akl 
1 19991 ) mounted on the Calar Alto 2.2-m telescope. BUSCA 


can obtain photometry in four different passbands simulta¬ 
neously, the incoming light being split by dichroics. Each 
of the four channels has a 4096 x 4096 pixel CCD, which 
were used with 2x2 binning to give a plate scale of 0.35" 
pixel -1 . The FOV of each channel depends on the filter in 
the light beam, being 5.8' in diameter for the Thuan-Gunn 
and Stromgren filters and 12' x 12' for the Cousins / filter. 
The telescope was autoguided and defocussed during both 
observing sequences. 

Unfortunately, the BUSCA controller requires the same 
exposure time to be used in all four channels. The exposure 
times were therefore chosen to avoid saturation in the r 
band, for which the count rate was highest. This meant 
that the redder channels yielded better-quality data than 
the bluer channels, especially the u band. 


2.4. Data reduction 


Data reduction was performed using standard methods. 
Bias and flat-field images on the sky were collected before 
and during twilight, respectively, and median-combined 
to generate master bias and flat-field frames. These were 
used to calibrate the science images. Light curves were 
then obtained usin g ape rture photometry algorithms from 
daophot (|Stetsonlll987h as implemented in the lDl0-based 
defot pipeline (see ISouthworth et al.ll2Ql3 and references 
therein), which uses subroutines from NASA’s astroli^I 
library. 

The sizes of the software apertures used for the aper¬ 
ture photometry are listed in Table[2j and were chosen to 


2 The acronym IDL stands for Interactive Data 

Language and is a trademark of ITT Visual 
Information Solutions. For further details see: 

http://www.ittvis.com/ProductServices/IDL.aspx 

3 The ASTROLIB subroutine library is distributed by NASA. 
For further details see: http://idlastro.gsfc.nasa.gov/. 
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Fig. 2. Light curves of transits of HAT-P-23b, observed 
with the Calar Alto 1.23-m telescope, compared with the 
best-fitting models given by jktebop. The date of each ob¬ 
servation is reported next to the corresponding light curve. 
The residuals from the fits are displayed at the base of the 
figure in the same order as the light curves. 
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be those that gave the lowest out-of-transit (OOT) scatter. 
We noticed that changes in the aperture size of both the 
target region and the sky annulus do not affect the over¬ 
all light-curve shape, but do cause a slight variation in the 
scatter of the data points. Once the aperture sizes were set, 
we extracted instrumental magnitudes for the target and 
possible comparison stars in the FOV. Pointing variations 
were corrected by cross-correlating each image against a 
reference frame. The light curves were then detrended by 
a second-order polynomial whilst optimising the weights of 
an ensemble of comparison stars. The choice of the compar¬ 
ison stars was carried out according to their brightness, by 
comparing the different OOT scatter and taking the com¬ 
bination that gave the lowest scatter. 
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Fig. 3. Light curves of transits of WASP-48 b, observed 
with the Cassini 1.52-m (“loiano”) and Calar Alto 1.23-m 
(“caha”) telescopes, compared with the best-fitting curves 
given by jktebop. The date of each observation is reported 
next to the corresponding light curve. The residuals from 
the fits are displayed at the base of the figure in the same 
order as the light curves. 


3. Light curve analysis 

We separately analyzed each of the _light curves observed 
using the JKTEBOF0 code (see ISouthworthll2008i and refer¬ 
ences therein), jktebop fits an observed light curve with a 
synthetic one, constructed according to the values of a set 
of initial parameters. Some of the input parameters are left 
free to vary until the best fit is reached. The main photo¬ 
metric parameters that we can measure with jktebop are 
the orbital inclination, i, the orbital period, P, the tran¬ 
sit midpoint, To, and the sum and ratio of the fractional 
radii of the star and planet, va + Od and k = v^/va- The 
fractional radii are defined as r\ = Pa/& and r b = Pb/&> 
where a is the orbital semi-major axis, and Ra and Pb are 
the absolute radii of the star and the planet, respectively. 

W e assumed circular orbits for both the plan etary sys¬ 
tems ([Enoch et al. 1 1201 ll: lO’Rourke et al.l 120141) . The val¬ 
ues of the planet-star mass ratio were fixed to the those 
obtained using the estimated masses from the discovery 
papers. The limb darkening effect on the light curves was 

4 The source code of jktebop is available at http:// 
www.astro.keele.ac.uk/jkt/codes/jktebop.html 
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Fig. 4. A planetary transit event of HAT-P-23b as ob¬ 
served with the BUSCA instrument mounted on the CAHA 
2.2-m telescope. The light curves in the Thuan-Gunn ugrz 
are shown from top to bottom and are compared to the 
best jktebop fits. The residuals are plotted at the base of 
the figure. 


taken into account by modeling it using a quadratic law, 
and we checked that the difference in the results obtained 
using a linear or a logarithmic law is negligible. 

The fits of the light curves were performed using the¬ 
oretical values of t he limb darkening coefficients ([Claretl 
120001: IClaret 1 120041) and by fitting the linear coefficients 
whilst fixing the quadratic coefficients but perturbing them 
during the error analysis simulations. The atmospheric pa¬ 
rameters of the stars HAT-P-23A and WASP-48 A, as¬ 
sumed for deriving the initial values of the limb-darkening 
coefficients, are listed in Table[3j as well as the weighted- 
mean values of the limb-darkening coefficient obtained from 
the fit of the light curves. 

Since the aper routine, which we used to perform aper¬ 
ture photometry on the calibrated science images, com¬ 
monly underestimates the error bars, we enlarged them for 
each dataset by multiplying the errorbar for each photome¬ 
tric point by the square-root of the xl obtained through an 
initial fit of the corresponding light curve. We then further 
inflat e d the errors using the /3 approach (e .g. iPont et all 
120061 : iGillon et al.l 12006c IWinn et all 120071) to take into 
account the presence of systematic effects and correlated 
noise. 


WASP-48 


1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 [ 1 1 1 I 
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-0.08-0.06-0.04-0.02 0.00 0.02 0.04 0.06 
Phase 

Fig. 5. As for Fig.0]but for WASP-48b. A Cousins / filter 
was used in the reddest arm of BUSCA, while a Stromgren 
u was used in the bluest one. 


Table 2. Stellar atmospheric parameters used to calculate 
the limb darkening coefficients, and weighted-mean values 
of the linear coefficients obtained from the fit of the light 
curves. 


Parameter 

HAT-P-23 A 

WASP-48 A 

T e g(K) 

5905 

6000 

log g (cm s -2 ) 

4.33 

4.03 

[f] 

0.1 

-0.1 

Umicro (km S ) 

2.0 

2.0 

LD coeff R 

0.317 ±0.031 

0.277 ± 0.075 

LD coeff I 

- 

0.101 ±0.067 

LD coeff u 

0.780 ±0.511 

1.09 ±0.33 

LD coeff g 

0.239 ± 0.090 

0.631 ±0.164 

LD coeff r 

0.221 ± 0.066 

0.229 ± 0.065 

LD coeff z 

0.179 ±0.139 

- 


To assign the uncertainties to each parameter obtained 
from the fitting process, we generated 10 000 simulations 
with a Monte Carlo algorithm and also used a residual- 
permutation algorithm (ISouthworthl 12008 b Since none of 
the two techniques systematically gives lower uncertainties, 
we took the larger of the 1-a values obtained using the 
two algorithms. The best-fitting model and the residuals 
for each of the 21 light curves are shown in Figures [2j [3j 0J 
and [5] 
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Table 3. Times of transit midpoi nt of H AT-P-23b and 
their residuals. References: (1) iBakos et al.l (2 0lf[) : (2) 
TRESCA; (3) CA 1.23-m, this work; (4) CA 2.2-m, this 
work. 


Time of minimum 
BJD(TDB) -2400000 

Epoch 

Residual 

(JD) 

Reference 

54852.265383 ± 0.00018 

0 

-0.00061 

(i) 

55434.453087 ± 0.00103 

480 

0.00331 

(2) 

55451.430268 ± 0.00112 

494 

0.00013 

(2) 

55736.460815 ± 0.00086 

729 

0.00321 

(2) 

55749.802945 ± 0.00072 

740 

0.00362 

(2) 

55753.439625 ± 0.00112 

743 

0.00166 

(2) 

55776.484166 ± 0.00085 

762 

0.00142 

(2) 

55783.762296 ± 0.00051 

768 

0.00226 

(2) 

55812.871847 ± 0.00114 

792 

0.00262 

(2) 

55838.341217 ± 0.00063 

813 

0.00145 

(2) 

56129.433378 ± 0.00060 

1053 

0.00172 

(2) 

56129.433908 ± 0.00068 

1053 

0.00225 

(2) 

56129.432399 ± 0.00041 

1053 

0.00057 

(3) 

56135.496648 ± 0.00157 

1058 

0.00263 

(2) 

56135.495288 ± 0.00032 

1058 

0.00158 

(3) 

56159.756359 ± 0.00085 

1078 

0.00202 

(2) 

56186.438739 ± 0.00090 

1100 

-0.00075 

(2) 

56460.547674 ± 0.00033 

1326 

0.00124 

(3) 

56477.528543 ± 0.00025 

1340 

-0.00137 

(3) 

56478.743946 ± 0.00091 

1341 

0.00074 

(2) 

56500.574007 ± 0.00281 

1359 

-0.00079 

(3) 

56505.425061 ± 0.00052 

1363 

-0.00101 

(3) 

56539.383814 ± 0.00168 

1391 

-0.00050 

(4) u 

56539.385859 ± 0.00036 

1391 

0.00019 

(4)3 

56539.385926 ± 0.00021 

1391 

-0.00029 

(4) r 

56539.385145 ± 0.00051 

1391 

-0.00226 

(4 )z 

56562.430097 ± 0.00122 

1410 

-0.00021 

(2) 

56562.432087 ± 0.00098 

1410 

-0.00015 

(2) 

56579.409838 ± 0.00114 

1424 

-0.00093 

(2) 


3.1. New orbital ephemeris 

We refined the orbital periods of both planets by using our 
new photometric data and transit timings available in the 
literature or at the ETE0 website. 

We made a linear fit to the transit midpoints versus 
cycle numbers in order to improve the ephemeris. All the 
transits considered in the linear fit for the two planetary 
systems are listed in Tables [3] and 0J The new values for 
the ephemeris found from the fits are: 

T 0 = BJD(TDB) 2 454 852.26599 (20) + 1.21288287 (17) E, 

and 

T 0 = BJD(TDB) 2 455 364.55241 (24) + 2.14363544 (58) E. 

for HAT-P-23 and WASP-48, respectively. The numbers in 
brackets are the errors relative to the last digit, and E 
represents the cycle number. 

If a TEP system, known to be composed of a parent 
star and a planet, hosts other planets, the gravitational 
interaction between the planetary components results in a 
periodic delay and advance in the times of transit of the 
known planet. We checked for a possible third component 
in the HAT-P-23 and WASP-48 systems by looking to see 
whether there was a periodic variation in the transit times 
of the known planets. 

The residuals from the linear fits are plotted in Figs. 
[6] and 0 as a function of cycle number and do not show 
any clear systematic deviation from the predicted transit 

5 The Exoplanet Transit Database (ETD) website can be 
found at http: //var2. astro. cz/ETD 


Table 4. Times of transit mid-p oint of WA SP -48 b and 
their residuals. References: (1) lEnoch et al. I ( 2011 b (2) 
TRESCA; (3) Loiano 1.52-m, this work; (4) CA 1.23-m, 
this work; (5) CA 2.2-m, this work. 


Time of minimum 
BJD(TDB)—2400000 

Epoch 

Residual 

(JD) 

Reference 

55364.55202 ± 0.00027 

0 

-0.00039 

(i) 

55696.80915 ± 0.00131 

155 

-0.00675 

(2) 

55696.81647 ± 0.00155 

155 

0.00056 

(2) 

55750.41680 ± 0.00151 

180 

0.01000 

(2) 

55782.56695 ± 0.00195 

195 

0.00562 

(2) 

55797.56545 ± 0.00160 

202 

-0.00131 

(2) 

55825.43480 ± 0.00152 

215 

0.00076 

(2) 

56033.37084 ± 0.00162 

312 

0.00416 

(2) 

56065.51850 ± 0.00187 

327 

-0.00270 

(2) 

56168.42213 ± 0.00067 

375 

0.00642 

(2) 

56391.34809 ± 0.00218 

479 

-0.00570 

(2) 

56393.49772 ± 0.00087 

480 

0.00029 

(2) 

56453.51930 ± 0.00181 

508 

0.00007 

(2) 

56468.53007 ± 0.00110 

515 

0.00540 

(2) 

56487.81919 ± 0.00060 

524 

0.00180 

(2) 

56511.40314 ± 0.00252 

535 

0.00576 

(2) 

56541.40894 ± 0.00103 

549 

0.00067 

(2) 

56541.41054 ± 0.00139 

549 

0.00227 

(2) 

55705.37900 ± 0.01013 

159 

-0.01145 

(3) 

55707.53160 ± 0.00041 

160 

-0.00248 

(3) 

55797.57056 ± 0.00063 

202 

0.00378 

(4) 

55797.56621 ± 0.00162 

202 

-0.00056 

(5) u 

55797.56614 ± 0.00063 

202 

-0.00063 

(5) 9 

55797.56473 ± 0.00053 

202 

-0.00204 

(5) r 

55797.56569 ± 0.00050 

202 

-0.00108 

(5) I 

56183.42469 ± 0.00062 

382 

0.00353 

(4) 

56498.53535 ± 0.00025 

529 

-0.00021 

(4) 

56811.50538 ± 0.00025 

675 

-0.00094 

(4) 


times. However, the quality of the fits, xl = 3.26 for HAT- 
P-23 and xl — 10.05 for WASP-48, indicates that a lin¬ 
ear ephemeris is not a good match to the observations in 
both the cases. Based on our e xperience with a similar s it¬ 
uation in previous st udies (e.g. ISouthworth et al.ll2012 a.b: 
iMancini et al.ll2013allcl ). we conservatively do not interpret 
the large xl values as sign of transit timing variations, but 
as an underestimation of the uncertainties in the various To 
measurements. 

3.2. Final photometric parameters 

The final photometric parameters for the HAT-P-23 sys¬ 
tem were obtained taking into account only the complete 
transit light curves (i.e. discarding the two partial ones). 
This choice was made after noticing that the x^, rela¬ 
tive to the ratio of the radii of the planet and the star, 
decreases significantly, from 4.37 to 0.99, when rejecting 
them. Concerning WASP-48, we discarded the incomplete 
light curve and those with very high scatter (i.e. the light 
curve observed on 23-08-2011 with the CA 1.23-m and with 
BUSCA in the u band). Even discarding three light curves, 
the Xu f° r k parameter remained rather high, 8.48. This 
discrepancy is discussed further in Section lCT 

The final values were obtained as a weighted mean of 
all the transits taken into account. The relative errors, ob¬ 
tained from the weighted mean, are rescaled by multiplying 
them for the relative x£- 

The parameters of the jktebop fits to each of our light 
curves are given in Tables [5] and [6] for HAT-P-23 and WASP- 
48, respectively. These tables also show the final photome- 
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Fig. 6. Plot of the residuals of the timing of mid-transit of HAT-P-23 versus a linear ephemeris. The points indicate 
literature results (blue), data obtained from the TRESCA catalogue (red), and our data (green). 
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Fig. 7. Plot of the residuals of the timing of mid-transit of WASP-48 versus a linear ephemeris. The points indicate the 
value from the discovery paper (blue), data obtained from the TRESCA catalogue (red), and our data (green). 


Table 5. Photometric properties of the HAT-P-23 system derived by fitting the light curves with jktebop. The final 
parameters are given in bold and are compared with those of the discovery paper. (*) Light curve not taken into account 
for the final result. 


Source 


r*A ± r h 

k 

i 

r*A 

Tb 

CA 1.23-m # 

1 * 

0.26447 ± 0.02107 

0.11350 ±0.00267 

89.72 ± 0.90 

0.23751 ± 0.01893 

0.02696 ± 0.00224 

CA 1.23-m # 

2 * 

0.28307 ±0.01135 

0.12220 ±0.00103 

82.01 ± 1.19 

0.25225 ±0.01002 

0.03083 ± 0.00138 

CA 1.23-m # 

3 

0.23870 ± 0.00861 

0.11343 ±0.00129 

88.48 ± 2.23 

0.21438 ± 0.00765 

0.02432 ± 0.00102 

CA 1.23-m # 

4 

0.24324 ± 0.00856 

0.11196 ±0.00180 

87.05 ± 2.44 

0.21875 ± 0.00745 

0.02450 ±0.00116 

CA 1.23-m # 

5 

0.25716 ± 0.00625 

0.11732 ±0.00113 

84.48 ± 0.92 

0.23016 ± 0.00539 

0.02700 ± 0.00086 

CA 1.23-m # 

6 

0.25561 ± 0.01734 

0.11175 ±0.00498 

86.56 ± 3.42 

0.22992 ± 0.01533 

0.02569 ± 0.00216 

CA 1.23-m # 

7 

0.25067 ± 0.01582 

0.12028 ±0.00201 

85.94 ± 3.22 

0.22376 ± 0.01382 

0.02691 ± 0.00205 

CA 2.2-m u 


0.28919 ± 0.03285 

0.12114 ±0.01224 

81.87 ±4.85 

0.25794 ± 0.02940 

0.03125 ± 0.00462 

CA 2.2-m g 


0.25315 ± 0.02248 

0.11856 ±0.00192 

84.18 ±3.71 

0.22632 ± 0.01972 

0.02683 ± 0.00278 

CA 2.2-m r 


0.24223 ± 0.00920 

0.11668 ±0.00121 

86.30 ± 2.50 

0.21692 ±0.00819 

0.02531 ±0.00124 

CA 2.2-m z 


0.23377 ± 0.00981 

0.11505 ±0.00207 

87.17 ±2.64 

0.20965 ± 0.00854 

0.02412 ±0.00131 

Final results 

0.24539 ± 0.00499 

0.11616 ±0.00081 

85.74 ±0.95 

0.21998 ±0.00436 

0.02541 ± 0.00065 

Bakos et al. f20111 

0.268 ±0.014 

0.1169 ±0.0012 

85.1 ± 1.5 

0.240 ± 0.014 

0.028 ±0.014 


Table 6. Photometric properties of the WASP-48 system derived by fitting the light curves with jktebop. The final 
parameters are given in bold and are compared with those of the discovery paper. (*) Light curve not taken into account 
for the final result. 


Source 

Cassini # 1 * 
Cassini # 2 
CA 1.23-m # 1 
CA 1.23-m # 2 
CA 1.23-m # 3 
CA 2.2-m u * 
CA 2.2-m g 
CA 2.2-m r 
CA 2.2-m I 
CA 1.23-m # 4 


r*A + r h 


k 


r a 


Tb 


0.23478 

0.22734 

0.18406 

0.23268 

0.23979 

0.28972 

0.20267 

0.22335 

0.21368 

0.23689 


± 0.02433 
± 0.01098 
± 0.00280 
± 0.01670 
± 0.00790 
± 0.04443 
±0.01931 
±0.02172 
± 0.01603 
±0.01140 


0.08349 

0.09028 

0.07810 

0.09426 

0.09648 

0.09503 

0.08247 

0.08900 

0.09392 

0.09613 


± 0.00233 
± 0.00081 
±0.00131 
±0.00201 
± 0.00046 
±0.01380 
± 0.00407 
± 0.00367 
± 0.00145 
±0.00126 


84.33 ± 2.50 

82.34 ± 0.96 
89.90 ± 0.65 
82.49 ±1.70 
81.56 ±0.67 
78.88 ± 3.66 
85.24 ± 3.25 
83.23 ± 2.25 
83.98 ± 1.63 
81.79 ±1.04 


0.21668 ±0. 
0.20852 ± 0. 
0.17072 ±0. 
0.21264 ±0. 
0.21870 ±0. 
0.26458 ± 0. 
0.18723 ±0. 
0.20509 ± 0. 
0.19533 ±0. 
0.21612 ±0. 


02273 

00995 

00270 

01496 

00715 

03772 

01750 

01958 

01449 

01020 


0.01809 ± 
0.01883 ± 
0.01333 ± 
0.02004 ± 
0.02110 ± 
0.02515 ± 
0.01544 ± 
0.01825 ± 
0.01835 ± 
0.02078 ± 


0.00184 

0.00102 

0.00022 

0.00179 

0.00074 

0.00678 

0.00193 

0.00230 

0.00156 

0.00123 


Final results 

0.23276 ±0.00578 

0.09584 ±0.00077 

81.99 ±0.54 

0.21272 ±0.00520 

0.02010 ±0.00059 

Enoch et al. (2011) 

0.259 ±0.013 

0.098 ±0.014 

80.09io.79 

0.2364 ±0.0125 

0.0227 ± 0.0272 
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Table 7. Spectroscopic parameters of the stars HAT-P- 
23 A and WASP-48 A. 


Source 

HAT-P-23 A 

Ref. 

WASP-48 A 

Ref. 

T eff (K) 

5885 ± 72 

1 

6000 ± 150 

4 

[f] 

0.13 ±0.08 

1 

-0.12 ±0.12 

4 

Ka (ms 

368.5 ± 17.6 

2 

136.0111;? 

4 

e 

0.0 fixed 

3 

0.0 fixed 

4 


Notes. (1) iTorres et al.l (EiFlTIl: (2) iBakos et alJ <f201lh : (3) 
IQ’Rourke et, all (I2014D : f4') lEnoch et, al. I f 201 lh . 


Table 8. Physical properties for the HAT-P-23 system. The 
results obtained in this work are reported and compared 
with those of the discovery paper. 



This work (final) 

Bakos et al. (2011) 

Ma (M©) 

1.104 ±0.043 ±0.018 

1.13 ±0.04 

Ra (Rq) 

1.089 ±0.027 ±0.006 

1.20 ±0.07 

log PA (cgs) 

4.407 ±0.018 ±0.002 

4.33 ± 0.06 

PA ( PQ) 

0.855 ±0.051 


M b (Mj U p) 

2.07 ±0.12 ±0.02 

2.060 r 0.1 1 1 

Rb (-Rjup) 

1.224 ±0.036 ±0.007 

1.368 ±0.090 

Pb (ms -2 ) 

34.3 ± 2.4 

27.5 ±3.2 

Pb (Pjup) 

1.057 ±0.097 ±0.006 

0.81 ±0.15 

r eq (K) 

1951 ± 30 

2056 ± 66 

e 

0.0706 ± 0.0040 ± 0.0004 

0.062 ± 0.004 

a (AU) 

0.02302 ± 0.00030 ± 0.00012 0.0232 ± 0.0002 

Age (Gyr) 

9 -1+3. 4+1.1 
z - 1 -4.5 -2.0 

4.0 ± 1.0 


trie parameters in bold font, and the results from the dis¬ 
covery papers for comparison. 


4. Physical properties of HAT-P-23 and WASP-48 


Following the methodology used bv ISouthworthl (l2010lb the 
main physical parameters of the TEP systems HAT-P-23 
and WASP-48 were found from the photometric param¬ 
eters deduced from the parameters available in the liter¬ 
ature (see Table0 for values and references), and by in¬ 
terpolating within the tabul ated predictions of five sets 
of th e oretical stellar models ([ Claret 2004: iDemaraue et al.l 
200 4 iPietrmferni et al.l 12004 : YandenBerg et al.l 120061: 
Potter et al. 2008). In brief, an initial estimate of the stel¬ 


lar mass was specified and the observed quantities were 
compared to the ones predicted by stellar models for this 
mass. The mass was then iteratively adjusted to find the 
best agreement between the observed and expected values. 

Since the radius versus mass relation varies according 
to the age of the star, the interpolation was performed for 
different ages of the system, starting from 0.1 Gyr until the 
end of the main-sequence lifetime of the star, in steps of 
lOMyr. The output set of physical parameters is the one 
that gives the best agreement between the predicted and 
the measured quantities. Separate sets of results were cal¬ 
culated using each of the five sets of theoretical model tab¬ 
ulations. 

The final values found are shown in Table[8] for HAT- 
P-23 and Table[9] for WASP-48. Most quantities have two 
errorbars, and in these cases the first is a statistical er¬ 
ror obtained by propagating the uncertainties on the input 
measurements, and the second is a systematic error which 
is the scatter of the results fro m each of the five different 
sets of theoretical models ("see ISouthworthi’[2QQ9 ). 


Table 9. Final results for the physical parameters of 
WASP-48 obtained in this work compared to those of the 
discovery paper. 



This work (final) 

Enoch et al. (2011) 

M a (M©) 

1.062 ±0.074 ±0.014 

1.19 

±0.05 

Ra (Rq) 

1.519 ±0.051 ±0.007 

1.75 

±0.09 

log PA (Cgs) 

4.101 ±0.023 ±0.002 

4.03 

±0.04 

PA (PQ) 

0.303 ± 0.022 

0.22 

±0.03 

M b (Mj up ) 

0.907 ± 0.085 ± 0.008 

0.98 

±0.09 

Rb (-Rjup) 

1.396 ±0.051 ±0.006 

1.67 

±0.10 

p b (ms -2 ) 

11.5 ± 1.1 

8.1 

± 1.1 

Pb (Pjup) 

0.312 ±0.037 ±0.001 

0.21 

±0.04 

Teq (K) 

1956 ± 54 

2030 ± 70 

0 

0.0406 ± 0.0036 ± 0.0002 


— 

a (AU) 

0.03320 ± 0.00077 ± 0.00015 0.03444 

± 0.00043 

Age (Gyr) 

n £+1.0 +0.6 

U ' U -4.9 -2.8 

7 . 9 +° 


4.1. Radius vs wavelength variation 

Photometric observations of planetary transit events 
through different filters allow us to measure the appar¬ 
ent radius of transiting planets in each passband, obtain¬ 
ing an insight of the compo s ition of their atmo s phere 
e.g [Sonthworth et al.l 1201214 I Copper wheat et all 20131: 


Nik olov et al.l 2013 : iMancini et al.l 2013(I ~ Narita et al.l 

20131 ; iChen et al. 1 120141) . We used our multi-band data to 
investigate possible variations of the radius of HAT-P-23 b 
and WASP-48 b in different optical passbands. 

We phased all the light curves collected with the same 
telescope and filter combination an d, following the strategy 
used by ([Sonthworth et alll2Q12bl) . we re-fitted them and 
the BUSCA ones with all parameters fixed to the final val¬ 
ues reported in Table[8]and[9j with the exception of k. This 
approach allows to remove sources of uncertainty common 
to all data sets, maximizing the accuracy of estimations of 
the planet/star radius ratio as a function of wavelength. 
Due to the very large uncertainty, the values of k measured 
in the u band were ignored. The results in the other bands 
are shown in Fig. [8] for HAT-P-23 b and Fig. [9] for WASP- 
48 b. The vertical bars represent the relative errors in the 
measurements and the horizontal bars show the full width 
at half-maximum (FWHM) transmission of the passbands 
used. Transmission curves of the adopted filters are shown 
in the bottom panel. 

For illustration, the results obtained for HAT-P-23 b 
are compared with three spectra calculated from one- 
dimensional model atmospheres by iFortnev et al. I (120101 ) 
for a Jupiter-mass planet with a surface gravity of g b = 
25 ms -2 , a base radius of 1.25 Rj up at lObar, and T eq = 
2000 K. The first model (red line) is run in an isothermal 
case taking into account chemical equilibrium and the pre¬ 
sence of strong absorbers, like TiO and VO. The second 
model (green line) is obtained omitting the presence of the 
strong absorbers. The last model (blue line) is obtained by 
artificially removing the strong absorbers, as in the previous 
case, but also increasing the H 2 /He Rayleigh scattering by a 
factor of 100. The values o f k for WASP-48 b wer e compared 
with another model from IFortnev et al. I (|2010l) similar to 
previous ones, but for a Jupiter-mass planet with a surface 
gravity of g b = 10ms -2 and without strong absorbers. 

The precision of the final results does not allow us to 
discriminate among the models. We can only note that, 
since the are no large variations of the planets’ radius at 
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Fig. 8. Variation of the planetary radius of HAT-P-23 b, 
in terms of planet/star radius ratio, with wavelength. The 
black points are from the transit observed with BUSCA, 
while the blue point is the weighted-mean results coming 
from the seven transits observed in Cousins R with the 
CA1.23-m telescope. The vertical bars represent the er¬ 
rors in the measurements and the horizontal bars show the 
FWHM transmission of the passbands used. The observa¬ 
tional points are compared with three synthetic spectra for 
a Jupiter planet with a surface gravity of g b = 25 ms -2 
and T eq = 2000 K. The synthetic spectra in green and blue 
do not include TiO and VO opacity, while the spectrum in 
red does, based on equilibrium chemistry. With respect to 
the model identified with the green line, the blue one has 
H 2 /He Rayleigh scattering increased by a factor of 100. An 
offset is applied to all three models to provide the best fit to 
our radius measurements. Transmission curves of the filters 
used are shown in the bottom panel. 


shorter wavelength within the experimental uncertainties, 
the atmospheres of the two planets should be quite trans¬ 
parent and not affected by large Rayleigh scattering. They 
are well-suited for further investigation at the optical wave¬ 
lengths with more precise instruments. 

4.2. Eclipsing Binary 

The eclipsing binary candidate NSVS 3071474 is located at 
sky position RA(J2000) = 19 24 03.821 and Dec(J2000) = 
55 27 33.33. It is sufficiently close to WASP-48 that it is 
present in the FOV of some of the transits that we mon¬ 
itored. We have extracted light curves of this object and 
confirm the nature of the system to be that of a contact 
eclipsing binary. 

We also obtained th e li ght curve observed by the 
Super WASP survey ([Pollacco et al.ll2006h , where it is called 
IS WASP J192403.81+552734.5. The phased and binned 
WASP data, and the light curves obtained from our frames 
are shown in Figs. [TUI and [TT1 respectively. 

By fitting both the WASP data, and our lightcurves 
with JKTEBOP we obtained a measurement of the period 
of the eclipsing binary of 0.458998509 ± 0.000000011 days. 


Fig. 9. Variation of the planetary radius of WASP-48 b, 
in terms of planet/star radius ratio, with wavelength. The 
black points are from the transit observed with BUSCA, 
while the blue points are the weighted-mean results com¬ 
ing from the two transits observed in Cousins R and the 
two in I with the CA1.23-m telescope. The vertical bars 
represent the errors in the measurements and the horizon¬ 
tal bars show the FWHM transmission of the passbands 
used. The observational points are compared with a syn¬ 
thetic spectrum for a Jupiter planet with a surface gravity 
of #b = 10 ms -2 and T eq = 2000 K, which does not in¬ 
clude TiO and VO opacity. An offset is applied to all three 
models to provide the best fit to our radius measurements. 
Transmission curves of the filters used are shown in the 
bottom panel. 

Fig. 10. Light curves of NSVS 3071474 (1SWASP 
J192403.81+552734.5) phased (black dots) and binned (red 
dots) observed by the WASP survey between 2008 and 
2010. 

NSVS 3071474 



0.0 0.2 0.4 0.6 0.8 1.0 

Phase 


5. Discussion and Conclusions 

In this paper we presented refined parameters for the two 
TEP systems HAT-P-23 and WASP-48, obtained from new 
photometric data of transit events. We also presented simul¬ 
taneous observations in different optical bands. Our princi¬ 
pal results are as follows. 

• We confirmed the mass value of HAT-P-23 b, but ob¬ 
tained a radius for it smaller by roughly 1.5cr. This 
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Fig. 11. Light curves of NSVS 3071474, the eclipsing bi¬ 
nary close to WASP-48. The curves in the two top panels 
were observed with the 1.52-m Cassini telescope while the 
two bottom ones were observed with the 1.23-m Calar Alto 
telescope. 
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TEP now occupies a more populated region of the mass- 
radius diagram and is no longer found to be one of 
the highly inflated transiting planets (see Fig. EE}. The 
mean density is therefore higher than what was found 
in the discovery paper. 

• We obtained improved estimates of the mass and radius 
of both the star and the planet of the WASP-48 plane¬ 
tary system. The values for the stellar and planetary 
radius are smaller than those found in the discovery pa¬ 
per by roughly 2.2cr and 2.4a, respectively. The masses 
are consistent with previous results. 

• A study of the planets’ radius variation as a function 
of optical wavelength, based on the data presented in 
this work, do not indicate any large variation for ei¬ 
ther planet, suggesting that their atmospheres are not 
affected by a large Rayleigh scattering. Further investi- 


Fig. 12. Mass versus radius diagram of the transiting plan¬ 
ets. The blue and red dots represent the values for HAT-P- 
23 and WASP-48 found in literature and obtained in this 
work respectively. 
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gations of the transmission spectra of these two planets 
are needed to validate this statement. 

• Finally, we also presented new light curves for the eclips¬ 
ing binary NSVS 3071474, refining the measurement of 
the orbital period to be 0.458998509 (11) days. 
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